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ABSTRACT 

C^H' The primordial abundance of deuterium (D/H) yields a measure of the density 

of baryons in the universe and is an important complement to determinations from 
cosmic microwave background (CMB) experiments. Indeed, the current small sample 
^ ' of high rcdshift D /H measurements from quasar absorption line studies are in excellent 

I agreement with CMB-derived values. Conversely, absorption line measurements of the 

Galactic D/H ratio in almost 50 stellar sightlines show a puzzlingly large scatter 
outside the Local Bubble which is difficult to explain simply by astration from the 
primordial value. The currently favoured explanation for these large variations is that 
D is differentially depleted relative to H in some parts of the local interstellar medium 
(ISM). Here, we test this scenario by studying the correlation between D/H and the 
, abundance of titanium, one of the most refractory elements readily observed in the 

• ISM. Previous work by Prochaska, Tripp & Howk (2005) found tentative evidence 

for a correlation between Ti/H and D/H based on seven sightlines. Here we almost 
triple the number of previous Ti measurements and include several sightlines with 
very high or low D/H that are critical for quantifying any correlations with D/H. 
f-*) I With our larger sample, we confirm a correlation between Ti/H and D/H at the 97% 

confidence level. However, the magnitude of this dependence is difficult to reconcile 
with a simple model of dust depletion for two reasons. First, contrary to what is 
v> , expected from local depletion rates, the gradient of the highly refractory Ti is much 

j_j ■ shallower than that observed for Fe and Si. Second, we do not observe the established 

' tight, steep correlation between [Ti/H] and the mean volume density of hydrogen. 

Therefore, whilst dust remains a plausible explanation for the local D/H variations, 
the abundances of at least some of the refractory elements do not provide unanimous 
support for this scenario. We also argue that the correlations of [Si/H], [Fe/H], and 
[Ti/H] with D/H are inconsistent with a simple infall model of low metallicity gas 
with approximately solar abundances as the dominant cause for D variations. 
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1 INTRODUCTION 

The abundances of the light elements (deuterium, helium 
and lithium) form one of the cornerstones of Big Bang nucle- 
osynthesis theory and can all be used as baryometers, yield- 
ing the density of baryons in the universe as a fraction of the 
closure density, fi^. Although the light element abundances 



* Based on observations made with ESO Telescopes at the 
Paranal Observatories under programme ID 076.C-0503(A) 
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are actually a direct probe of the baryon-to-photon ratio, 
rj, this can be translated to Qt using the (now well-known) 
Cosmic Microwave Background (CMB) temperature. Deu- 
terium is by far the most sensitive of the light elements to 77, 
and thus the baryon density, and determining the primordial 
value of D/H (the deuterium to hydrogen abundance) has 
been the focus of intense research efforts over the last three 
decades (e.g. York & Rogerson 1976; Vidal-Madjar et al. 
1977; Laurent, Vidal-Madjar & York 1979; Pettini & Bowen 
2001; O'Meara et al. 2001; Hebrard et al. 2002; Moos et al. 
2002; Sonneborn et al. 2002; Kirkman et al. 2003; Sembach 
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et al. 2004; Wood et al. 2004; Oliveira et al. 2006). These 
efforts have made a two-pronged attack, by measuring D /H 
abundances both in the Milky Way and at high redshifts. 
Although D is not produced in any significant quantities af- 
ter Big Bang nucleosynthesis, it is easily astrated in stars, 
so measurements at high redshifts could potentially yield 
the primordial abundance more readily than observations 
in the local universe. High redshift measurements can be 
achieved with observations of high HI column density ab- 
sorption systems in the lines of sight towards QSOs, such as 
the damped Lyman alpha systems (DLAs) and Lyman limit 
systems (LLSs). However, since the DI Lyman a line is only 
separated from that of HI by 80 km/s, only those rare 
absorbers with very simple velocity structure yield reliable 
results. Thus, despite a decade of high redshift work, only 
six robust measurements of D/H at high redshift currently 
exist (O'Meara et al. 2001; Pettini & Bowen 2001; Kirkman 
et al. 2003; O'Meara et al. 2006). Nonetheless, these mea- 
surements paint a relatively coherent picture and yield an Qt 
(inferred from a mean D/H = (2.8±0.4) xl0"^f|that is in 
solid agreement (i.e. within the 1 a errors) with that deter- 
mined from CMB experiments such as WMAP (see Pettini 
2006 for a review). 

Measurements of D/H in the Milky Way complement 
high redshift observations by providing an order of magni- 
tude more sightlines in which abundances can be measured 
(see Linsky et al. 2006 for the most recent compilation). 
Although some astration has probably occurred, the high- 
est Galactic D/H ratios set a lower limit to the primordial 
abundance. Moreover, modelling the amount of astration in- 
ferred from sub-primordial abundances provides clues to the 
rate at which material has been processed in stars - a vital 
ingredient in chemical evolution models. For these reasons, 
one of the principle drivers of the Far-Ultraviolet Spectro- 
scopic Explorer (FUSE) was to compile a large database of 
Galactic D/H measurements: 'The main science goal of the 
FUSE mission has been to obtain accurate measurements 
of (D/H) gas for many sightlines in the Milky Way Galaxy 
and beyond in order to measure (D/H)prim and to obtain 
constraints on Galactic chemical evolution' (Linsky et al. 
2006). However, despite a massive investment of telescope 
resources from both FUSE and other satellite facilities, and 
measurements of D/H in almost 50 lines of sight, the picture 
of Galactic deuterium abundances remains puzzling. On the 
one hand, inside the Local Bubble (< 100 pc from the Sun) 
the D/H value seems roughly constant at 1.5 xlO~^, ap- 
proximately 60% of the primordial value. However, beyond 
this bound, there is a large scatter in the D/H ratios which 
can not be easily explained by astration (Jenkins et al. 1999; 
Wood et al. 2004; Draine, 2006). As pointed out by Draine 
(2006), there is an observed factor of five range in D/H for 
clouds within 100 pc. It is difficult to reconcile such dis- 
parate abundances, on such small spatial scales, within a 
chemical evolution model, particularly since mixing should 
be efficient on these scales (e.g. Romano et al. 2006). 

At the present time, the most favoured explanation for 
the variation in Galactic D /H ratios is the depletion of deu- 



^ The error on the D/H abundance from QSO absorption lines is 
estimated from a jackknife analysis of the weighted means of six 
robust measurements. 



terium from the gas phase onto grains. Deuterium depletion 
was first suggested over 20 years ago by Jura (1982) and 
then explored theoretically by Tielens (1983). The most re- 
cent theoretical work has been done by Draine (2004, 2006) 
who suggested that deuterated polycyclic aromatic hydro- 
carbons (PAHs) in cold dust could be sufficient sinks for 
deuterium to explain the scatter in Galactic measurements. 
Empirical support for D depletion comes from the compo- 
sition of interstellar grains collected in the upper terrestrial 
atmosphere which can be enriched by D relative to the usual 
meteoritic composition (Keller, Messenger & Bradley 2000). 
Whilst the quantity of deuterium was not sufficient to ex- 
plain the Galactic ISM abundances, the presence of dust 
phase D does act as a 'proof of concept'. 

One empirical test of the depletion scenario is to look 
for correlations between D/H and the abundances of other 
refractory elements. This test was first done by Prochaska, 
Tripp & Howk (2005) using far-blue optical lines of Ti II 
measured in Keck/HIRES spectra for seven Galactic sight- 
lines. Prochaska et al. (2005) noted a possible correlation at 
the 95% level, although this was strongly driven by a single 
high D/H sightline. Moreover, Prochaska et al. note that the 
literature limit for a second high D /H sightline (7^ Vel) does 
not agree with the overall trend in their data and that there 
is a large scatter in Ti/H at low D/H. More recently, Linsky 
et al. (2006) have presented corroborating evidence for de- 
pletion correlations based on Si and Fe; both elements show 
significant abundance trends with D/H. Linsky et al. also 
tentatively suggest that a possible correlation between D /H 
and the gas temperature derived from H2 may also support 
the dust scenario. Savage et al. (2007) have also supported 
the dust depletion model, since their high value (2.2xl0~^) 
of D/H in the warm neutral medium of the Galactic halo 
is consistent with expulsion of material from the disk and 
associated ablation of dust particles. 



In this paper, we extend previous analyses of refractory 
elements by significantly enlarging the sample of Ti mea- 
surements in the local ISM. Due to its highly refractory na- 
ture (e.g. Savage & Sembach 1996; Jenkins 2004), Ti should 
provide one of the most sensitive measurements of depletion 
and therefore correlate most steeply with D/H. A further 
advantage of Ti II is its ionization potential (IP) of 13.6 eV 
which means that ionization corrections will be less than for 
other commonly detected species such as Fe II (IP=16.2 eV) 
and Si II (IP=16.3 eV). 



2 OBSERVATIONS AND DATA ANALYSIS 

We chose targets with measured D/H and whose log N(HI) 
> 19.7. This criterion ensures that Ti II lines will have de- 
tectable equivalent widths (EWs) and that ionization cor- 
rections will be a negligible concern. 

2.1 Observations and Data Reduction 

Our final sample of Ti II lines in the local Galactic disk 
comprises data obtained from the following sources: the pi- 
lot survey of Prochaska et al. (2005) (7 stars), VLT/UVES 
data obtained from the archive (2 stars), VLT/UVES data 
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Table 1. Journal of Observations 



Star Mag (V) RA (J2000) Doc (J2000) Distance (pc) Instrument Run 



5 Ori 


2.23 


05 


32 


00, 


.4 


—00 17 57 


281+65 


HIRES 


September 2004 


i Ori 


2.77 


05 


35 


26 


.0 


—05 54 35.6 


407+127 


HIRES 


September 2004 


e Ori 


1.70 


05 


36 


12 


.8 


— Ui iZ Uu.y 


A 1 1 c: A 




beptember 2UU4 


Col 


5.17 


05 


45 


59 


.9 


— 62 io 26.2 


/I nn+ 100 


UVES 


Archive 


HD 411d1 


6.76 


06 


05 


52 


.5 


+48 14 57.4 


1253 


HIRES 


April 5 2006 


HD 53975 


6.47 


07 


06 


36, 


.0 


-12 23 38.2 


1318 


HIRES 


April 5 2006 


CPup 


2.25 


08 


03 


35, 


.1 


-40 00 11.3 


429+94 


UVES 


Period 76 


7^ Vel 


1.78 


08 


09 


32, 


.0 


-47 20 11.7 


258+35 


UVES 


Period 76 


WD 1034+001 


13.22 


10 


37 


04, 


.0 


-00 08 20 




HIRES 


June 2006 


e Car 


2.76 


10 


42 


57, 


.4 


-64 23 40.0 


135+9 


UVES 


Period 76 


a Cru 


1.33 


12 


26 


35, 


.9 


-63 05 56.7 


98+6 


UVES 


Period 76 


13 Cen 


0.61 


14 


03 


49, 


.4 


-60 22 22.9 


161+16 


UVES 


Archive 


BD+39°3226 


10.18 


17 


46 


31, 


.9 


+39 19 09.07 


290^1^° 


HIRES 


April 5 2006 


HD 191877 


6.26 


20 


11 


21, 


.0 


+21 52 29.8 


2200+550 


HIRES 


September 2004 


HD 195965 


6.98 


20 


32 


25, 


.6 


+48 12 59.3 


794+200 


HIRES 


September 2004 


HD 206773 


6.9 


21 


42 


24, 


.2 


+57 44 10 


500l?oo 


HIRES 


June 2006 


BD+28°4211 


10.51 


21 


51 


11, 


.0 


+28 51 50.4 


104+18 


HIRES 


October 2004 


HD 209339 


6.69 


22 


00 


39, 


.3 


+62 29 16.0 


850 


HIRES 


June 2006 


Feige 110 


11.83 


23 


19 


58, 


.4 


-05 09 56.2 


179+265 


HIRES 


September 2004 



obtained through a Period 76 ESO allocation (4 stars) and 
new HIRES observations obtained in 2006 (6 stars). For 
completeness, we list all targets from our Ti II survey (which 
include those already published in Prochaska et al. 2005) in 
Table [T] but only describe the data acquisition of previously 
unpublished data in the following sections. The distances 
to the stars are taken from the compilation of Linsky et 
al. (2006), supplemented by values from Pan et al. (2004), 
Oliveira & Hebrard (2006) and Hebrard et al. (2005a). 



2.1.1 HIRES Data 

The new HIRES data were obtained, primarily during 
twilight, on the nights of 06 April 2006, 02 June 2006, 
03 June 2006, and 25 December 2006. All of the data 
were acquired with the 0.4" wide, E3 decker which af- 
fords a FWHM resolution of « 3km s~^. We reduced 
these data using the HIRedux pipeline (v. 2.1; Bern- 
stein et al., in preparation) developed using the IDL 
software package and distributed within the XIDL pack- 
age (http://www.ucolick.org/~xavier/IDL/index.html). In 
brief, the pipeline bias subtracts and flatfields each science 
and calibration image. It traces the order curvature across 
each chip (limited to the blue CCD for these observations) 
and derives a 2D wavelength solution from the ThAr images. 

For these very bright targets, the biggest background is 
scattered light. This component is fit with a 2D surface by 
measuring the signal between echelle orders. The fit is in- 
terpolated across the orders and subtracted from the image. 
The RMS error in this signal is typically less than 5 elec- 
trons. We then boxcar extracted the stars (the S/N is too 
large for optimal extraction to be advantageous) to derive a 
ID spectrum and co-added the individual exposures, weight- 
ing by the median S/N. 



2.1.2 Archival UVES Data 

Data for /i Col and (3 Cen were downloaded from the ESO 
archive. A considerable amount of data were available but 
we selected a subset whose settings offered the highest res- 
olution and maximum exposure times at the wavelengths of 
interest. For j3 Cen we used only 1x1 binned data with 4 
central wavelengths: 346 nm (slit width = 0.5 arcsec), 390 
nm (slit width = 0.44 arcsec), 437 nm (slit width = 0.5 arc- 
sec) and 564 nm (slit width = 0.3 arcsec) yielding a FWHM 
resolution ranging from R ^ 70, 000 to 90,000 depending 
on the slit and the CCD. For ^ Col, the data were more 
limited and often had sub-optimal instrument set-ups (such 
as very wide slits). We therefore used only 1 central wave- 
length setting, 390 nm, with 2x2 binning and a 1 arcsec slit 
yielding a FWHM resolution R 45, 000. Although a set 
of MIDAS scripts exists to form a UVES pipeline, we found 
that the pipeline reduction had considerable trouble deal- 
ing with very high S/N data (a well-known problem). For 
consistency, we therefore decided to apply the HIRedux IDL 
reduction pipeline to the UVES data. This process worked 
well in general, although the final extracted spectra had con- 
siderable structure in the continuum (perhaps owing to poor 
flat-fielding), particularly in the blue orders. We attempted 
to fit out this structure during the continuum fitting process. 
This process was judged to be quite successful, based on the 
consistency of EWs of a given interstellar species between 
different orders, settings and transitions (see next section). 



2.1.3 UVES Period 76 Data 

For 4 southern stars not available in the archive, we received 
an ESO allocation in Period 76 (October 2005 - March 
2006). In all cases we used the same combination of set- 
tings: 1x1 binning, 0.5 arcsec slit and 3 central wavelengths, 
346 nm, 390 nm and 580 nm yielding a typical resolution 
R ~ 80,000. As for the archival UVES data, we used the 
adapted version of the HIRES IDL pipeline and fitted out 
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Figure 1. Normalized sections of spectra around the transitions 
of Ti II A 3384 and Ca II A 3934 for 4 of our program stars 
(see also Figure [2] and [Sjl . Note the different y-axes from panel to 
panel. 



Figure 2. Normalized sections of spectra around the transitions 
of Ti II A 3384 and Ca II A 3934 for 4 of our program stars (see 
also Figures [T] and [Sjl . Note the different y-axes from panel to 
panel. 



the structure in the blue orders along with the underlying 
continuum. 

3 RESULTS 
3.1 Analysis 

To determine gas phase column densities we use the VPFIlO 
software to decompose the absorption into individual com- 
ponents. The total column density is then calculated from 
the sum of these components. Profiles of Ti II A 3384, plus 
Ca II A3934 are shown in Figures [T] [2] and |3] In many lines 
of sight, Ti II A 3242 is also detected, in which case we si- 
multaneously fit the two transitions. The very small EWs 
(typically 10-20mA) of the Ti II hues coupled with the b- 
values (which indicate that the lines are resolved) show that 
it is unlikely that saturation affects these transitions. This is 
confirmed by an optical depth analysis of the two Ti II lines 
which yields consistent results for both. In the UVES data 
and a small subset of the HIRES data, Ca II A 3969 is also 
covered so that, when detected, it is simultaneously fit with 
Ca II A 3934. In the case of WD1034-h001, the Ti II and 

^ http:/ /www. ast.cam.ac.uk/~rfc/vpfit. html 



Ca II lines are barely detected, so the fit is very uncertain. 
In this case, we therefore additionally estimate the column 
density from the apparent optical depth of the absorption 
lines. For jj, Col, the Ti II A3384 line is detected at < 3a 
significance, so we quote a limit. We determine the errors 
on the lines in two different ways. First, we consider the for- 
mal error on the total column density returned by VPFnjf]. 
Next, we separately fit the different transitions of the same 
species, as well as fitting lines of a given species that are 
covered by multiple orders, or multiple settings. This is par- 
ticularly useful for the assessment of systematic errors such 
as the continuum fit, a particular concern for the blue UVES 
data. However, we found that the separate fits exhibited ex- 
cellent agreement, usually within 0.02 dex. Due to the very 
high S/N of the spectra, the range of column densities deter- 
mined from individual line fits usually dominated the formal 
error from VPFIT, so that the former was adopted as the 
final quoted uncertainty in Tabled We exclude two of the 
sightlines from our analysis of depletion patterns. First, HD 
206773 has poorly determined D /H due to severe saturation 
of the D I lines (Hebrard et al. 2005a). Second, although 



The errors on individual component column densities are also 
returned, but these are considerably more unreliable. 
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Table 2. Column Density Data 



Star 


N(H) 


N(D) 


N{OI) 


N{TiIl) 


N{CalI) 


N(Fell) 


N{SiII) 


Refs 


S Ori 


20.19±0.03 


15.06±0.12 


16.67±0.05 


11.15±0.04 




14.08±0.03 




1, 2, 3, 4 


L Ori 


20.18««? 


1 5 •^n+''-''"' 

i0.ciU_g 05 


16.76±0.09 


11.31±0.03 




14 20+'^'^° 




1, 3, 4, 5 


e Ori 


20.45±0.08 


15.25±0.05 


16.98±0.05 


11.40±0.03 




14.20±0.10 


15.02lO;?« 


1, 3, 4, 5 


Col 


19.86±0.02 


14 7+0-3 
' -0.1 




<11.4 


12.13±0.02 


14.13±0.02 


15.10±0.02 


4, 


6 


HD 41161 


21.08±0.09 


16.41±0.05 


18.04±0.05 


12.28±0.04 


12.54±0.01 


14.98±0.06 




7 




HD 53975 


21.14±0.06 


16.15±0.07 


17.87±0.08 


12.13±0.04 


12.56±0.01 


14.84±0.05 




7 




CPup 


19.96±0.03 


15.11±0.06 




11.07±0.04 


11.83±0.02 


14 13+0'17 




4, 


8 


7^ Vel 


19.71±0.03 


15.05±0.03 




li.lOiO.Oz 


"11 r--\ In no 

li.oiiO.Oz 


1 o n cr 1 n 1 1 

13.95±0.il 


1 cr no 1 n n ^1 

15.03±0.04 


4, 


8 


WD 1034+001 


20.07±0.07 


15.40±0.07 


16.60±0.10 


11.1±0.2 


11.35±0.15 


14.10±0.10 




4, 


9 


e Car 


20.28±0.10 


14 98+°-^* 
ii.a(5_0 21 


16.42±0.20 


11.31±0.01 


11.65±0.01 


14.19±0.03 


14.38l0;g5 


4, 


10 


a Cru 


19.85±0.10 


14.95±0.05 




11.12±0.02 


11.23±0.05 


14.00±0.10 


4, 


6 


P Cen 


19.63±0.10 


14.7±0.2 




11.05±0.04 


11.12±0.04 


, o 00+0.05 
io.yz_o 04 


14.59±0.01 


4, 


6, 11 


BD+-39°3226 


20.08±0.09 


15.15±0.05 


16.40±0.10 


11.49±0.08 


11.67±0.01 


14.15±0.07 


14.80±0.20 


4, 


9 


HD 191877 


21.10±0.07 


15 05+0-11 

iO.MO_Q og 


17 54+0.20 
l'-O+_0.12 


12.24±0.02 




14.95±0.02 




3, 


4, 12, 13 


HD 195965 


21.01± 0.05 


15.89±0.07 


.„ „_+0.04 
'-0.06 


12.02±0.02 




14.81±0.01 




3, 


4, 12, 13 


HD 206773t 


21.10±0.05 


17.0±0.5 


17.85±0.05 


12.16±0.02 


12.65±0.02 






14 


BD +28°4211 


19.85±0.02 


14.99±0.03 


16.22±0.10 


11.08±0.08 




14.10±0.10 




3, 


4, 15 


HD 209339tt 


~21.5 






12.14±0.03 


12.68±0.03 






16 




Feige 110 


20.14±0.07 


15.47±0.06 


17.06±0.15 


11.59±0.03 








3, 


17, 18 



Upper limits are 3cr. t Included here for completeness, but not used in analysis due to saturation of DI lines (Hebrard et al. 2005a). tt 
Included here for completeness, FUSE deuterium data still under analysis (C. Oliveira, private communication). Refs: 1 - Meyer et al. 

(1998); 2 - Jenkins et al. (1999); 3 - Prochaska et al. (2005); 4 - Linsky et al. (2006); 5 - Laurent et al. (1979); 6 - York & Rogerson 
(1976); 7 - Oliveira & Hebrard (2006); 8 - Sonneborn et al (2000); 9 - Oliveira et al. (2006); 10 - Allen et al (1992); 11 - van Steenberg 
& ShuU (1988); 12 - Hoopes et al. (2003); 13 - Liszt (2006); 14 - Hebrard et al. (2005a); 15 - Sonneborn et al. (2002); 16 - C. Oliveira, 
private communication; 17 - Friedman ct al. (2002); 18 - Hebrard et al. (2005b) 



HD 209339 has FUSE data which will eventually permit an 
analysis of the D/H in this sightline, the data have not yet 
been analysed. We tabulate the column densities of Ti II and 
Ca II for completeness, but do not consider these sightlines 
further in our analysis. 

For metal line species, we quote in Table [2] the col- 
umn densities of the single atomic or ionic species that is 
measured in our (or literature) data. However, Liszt (2006) 
has shown that the molecular contribution to D and H can 
be important and that ignoring its contribution can lead 
to an over-estimate of the D/H ratio (see also Lacour et 
al. 2005). We therefore use total column densities for hy- 
drogen and deuterium, accounting for both molecular and 
atomic gas, i.e. N(H)=N(HI)H-2N(H2)-f N(HD) and N(D) = 
N(DI)-|-N(HD). For most sightlines, the molecular contri- 
bution is orders of magnitude less than the atomic column 
density such that N(H) - N(H I) and N(D) ~ N(D I). How- 
ever, a significant effect is present for HD 41161, HD 53975, 
HD 191877, HD 195965, HD 206773 and HD 209339 (al- 
though these latter two are not used in our main analysis). 
The largest effect is seen towards HD 41161 where D I/H I = 
25.1 parts per million (ppm) and D/H = 21.4 ppm (Oliveira 
& Hebrard 2006), a difference of - 15%. 

3.2 Correlations with D/H 

In Figure|4]we present the first main result of this anal- 
ysis. The abundance of Ti, on both linear (in units of parts 
per billion, ppb) and logarithmic (relative to solar) scales is 
plotted versus D/H in units of ppm. For logarithmic abun- 
dances, we use the standard notation [X/H] = log (X/H) 



— log (X/H)q, which is analogous to the notation D(X) 
used by Linsky et al. (2006). For solar abundances, we use 
the compilation of Lodders (2003), except for oxygen, for 
which we adopt the value from Holweger (2001) where a 
full non-LTE treatment was adopted. These values are log 
[N(X)q/N(H)0] = -4.46, -4.53, -7.08, -3.26 for Si, Fe, Ti 
and O respectively. For convenience, in Table [3] we show 
the relative abundances used for Figures |4] and [5] We in- 
clude data from Prochaska et al. (2005) and perform a least- 
squares fit of the form y — mx -|- c to the final Ti sample 
which consists of 16 detections. The best fit straight line for 
[Ti/H] has values m = 0.016 ± 0.007 and c = -1.95 ± 0.09 
where errors were determined by randomly re-generating 
the [Ti/H] vs. D/H distribution 1000 times based on the 
quoted 1 a errors. Spearman's rank correlation statistic is 
0.55, yielding a rejection of the null hypothesis (that there is 
no correlation) at the 97% level. In addition to the enlarge- 
ment of the sample size, one of the important contributions 
of this work is the inclusion of 3 new high D/H data points, 
as well as the lowest D/H in the current sample. The correla- 
tion of Prochaska et al. (2005) was pivotal upon a single high 
datum (Feige 110); removal of this one point reduces the sig- 
nificance of the correlation from 95% to 86%. Moreover, for 
one other high D/H sightline in the literature, Prochaska 
et al. (2005) pointed out that the upper limit for 7^ Vel 
(D/H=21.9 ppm) of [Ti/H]s; -1.82 reported by Welsh et al. 
(1997) was potentially inconsistent with the depletion trend. 
With a larger sample, and more high D/H points, we can 
now make a more robust assessment of the D/H correlation. 
First, we find that with higher S/N data, the measured value 
of [Ti/H]=— 1.53 for 7^^ Vel is in actually excellent agreement 
with the high D/H point (Feige 110) presented in Prochaska 
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Table 3. Abundances 
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^^•^-4.4 


-1.89 ±0.21 
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-1.47 + 0.08 
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Figure 3. Normalized sections of spectra around the transitions 
of Ti II A 3384 and Ca II A 3934 for 4 of our program stars 
(see also Figure [T] and [2]| . Note the different y-axes from panel to 
panel. 



et al. Moreover, the highest D/H sightline (HD41161) in our 
sample is also in good agreement with the general trend of 
increasing Ti/H with D/H. However, the third high D/H 
sightline, WD1034+001 exhibits a low Ti/H for its D/H, al- 
though this point has the largest error bars on its measured 
N(Ti II) column density due to the low significance of the 
detection (see Figure |3] and Tables [2] and |3]) . 

Given the discrepancy noted above between the litera- 
ture upper limit and our measured value for [Ti/H] towards 
7 Vel, it is interesting to compare other sightlines for which 
literature values also exist. Welsh et al. (1997) report a limit 
on the column density of Ti II towards 6 Car that is half of 
our measured value and the discrepancy towards Pup is 
even higher: 0.45 dex, almost a factor of three. However, 
the Ca II abundances in the sightlines common to this work 
and Welsh et al. (1997) are in good agreement within the 
errors, with differences ^ 0.08 dex. Hunter et al. (2006) have 
also previously noted consistency in N(Ca II) and discrepan- 
cies in N(Ti II) between high S/N UVES data and previous 
works, whose S/N ratios are typically < 50. Hunter et al. 
attributed this to the weakness of the Ti II lines, relative 
to Ca II. Confirmation of this effect here highlights the im- 
portance of very high S /N ratios and spectral resolution to 
make these measurements accurately. 



4 DISCUSSION 
4.1 Dust Depletion 

Previous authors have interpreted the correlation between 
the gas-phase abundance of refractory elements and (D /H) 
as evidence that differential depletion leads to low values for 
each (Prochaska et al. 2005; Linsky et al. 2006). Although 
the trend traced by our enlarged sample of Ti/H values sup- 
ports this interpretation qualitatively, a closer examination 
of the results leads to a less clear cut picture. Focus first on 
the abundance trends from different elements. In Figure |3] 
we overlay the fits for depletion of Fe and Si from Linsky 
et al. (2006) with our own Ti results. If dust depletion is 
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Figure 4. Ti abundance as a function of D/If. Results from this 
work are sliow by stars; filled circles are data points taken from 
Prochaska et al. (2005). The solid line shows the least squares best 
fit to all 16 detections. The shaded region shows the constraints on 
the primordial value of D/H determined from quasar absorption 
line systems. In the top panel, we show linear-linear axes for Ti/H 
(in parts per billion, ppb) and D/H (in parts per million, ppm), 
the best fit line y = mx + c has m=0.06±0.03 and c=0.85±0.32. 
In the lower panel, we us a logarithmic y-axis to facilitate com- 
parison with other elements, the best fit line y = mx + c has 
m=0.016±0.006 and c = -1.95 ± 0.08. The fits to Fe and Si 
abundances from Linsky et al. (2006) are shown with dashed and 
dotted lines respectively. 



the cause of the correlation of abundance with D/H then 
we expect a steeper gradient for more refractory elements. 
In the case of Ti, Fe and Si, the former is the most re- 
fractory and the latter the least sensitive to depletion (e.g. 
Savage & Sembach 1996). The relative gradients in Figure |4] 
are therefore inconsistent with our simple predictions from 
a 'standard' depletion scenario, since we would expect Ti to 
be considerably steeper than Fe or Si. Indeed, the gradient 
seen in Ti is only marginally steeper than that found for 
[O/H] which is only very mildly refractory. Therefore, al- 
though the existence of trends between refractory elements 
and D/H support a picture of dust depletion, quantitatively 
the gradients are inconsistent with this scenario. Steigman, 
Romano & Tosi (2007) have also questioned the dust de- 
pletion scenario due to their failure to find a trend between 
line of sight reddening and D/H, as one might have expected 
were dust the cause of lower D/H values. 

In section 13.11 we argued that saturation of the Ti II 



lines is unlikely to be a serious concern. However, some of 
the Fe II and Si II transitions have considerably higher opti- 
cal depths (e.g., Redfield & Linsky 2002). There is an anti- 
correlation between [Fe/H] and N(Fe II) in the compilation 
of Linsky et al. (2006). If saturation systematically under- 
estimates N(Fe II) for the low [Fe/H] data points, the true 
slope of the correlation with D/H will be flatter than ob- 
served. However, this effect is unlikely to be large enough 
to yield a true gradient of [Fe/H] versus D/H that is flat- 
ter than the [Ti/H] versus D/H correlation. Moreover, there 
is a mild correlation between [Si/H] and N(Si II), so that 
any correction for saturation would make the dotted line in 
Figure |4] steeper, thus accentuating the difference with the 
[Ti/H] correlation. It is also possible that the D lines may 
suffer from saturation, and this may be particularly relevent 
for the high N(HI) lines of sight studied here. However, it 
is difficult to apply horizontal shifts to the D/H values in 
our sample that would make the correlation with [Ti/H] 
steeper. It therefore seems unlikely that saturation effects 
can explain the relative slopes in Figure |4] 

We can further investigate depletion trends by consid- 
ering correlations with mean hydrogen densities. It is well 
known that depletion in the Galactic ISM correlates with 
< n(H) >, the mean volume density of hydrogen (e.g. Sav- 
age & Bohlin 1979; PhiUips, Gondhalekar & Pettini 1982; 
Spitzer 1985; Gondhalekar 1985; Jenkins et al. 2004) where 



< n{H) >-- 



N{HI) + 2N{H2) 
distance 



(1) 



We calculate < n(H) > for the 17 sight lines in our Ti II 
sample and supplement them with five more sight lines in 
Linsky et al. (2006) which have H2 measurements (Lan 23, 
TDl 32709, PG0038-f 199, LSS 1274, and HD 90087) see Ta- 
bleland Figure [5] Jenkins (2004) parametrized ISM abun- 
dances in terms of a depletion factor, , which in turn cor- 
relates with < n(H) >. We performed a fit to the correla- 
tion between F-t and < n(H) > and using the parameters 
in Jenkins (2004), determined the corresponding relation- 
ship between < n(H) > and abundance for the 144 Galactic 
stars in his sample. For Ti, Si and Fe we determine; 

[Ti/H] = -2.226(0.4 < n(H) > +0.8) - 0.844 - 0.01 (2) 

[Si/H] = -1.076(0.4 < n(H) > +0.8) - 0.223 (3) 

[Fe/H] = -1.198(0.4 < n(H) > +0.8) - 0.95 + 0.02 (4) 

The final term in the equations for [Ti/H] and [Fe/H] ac- 
count for the slightly different solar abundance scale adopted 
by Jenkins (2004). This relationship between < n(H) > and 
ISM abundances is over-plotted with our data in Figure (5) 
The correlation between [Ti/H] and depletion factor in Jenk- 
ins (2004) is very steep (the steepest of the 15 elements he 
studies), and very tight (reduced x^=0.95). The best least 
squares fit to our Ti data is shown by the dashed line in the 
top left panel of Figure [S] The slope is clearly much flatter 
than that reported Jenkins (2004) and flatter than both the 
least squares fits for Fe and Si. Moreover, the correlation be- 
tween [Ti/H] and < n(H) > is not statistically significant: 
the Spearman rank correlation coefficient rules out the null 
hypothesis at < 2a (83%) significance. It is also interesting 
to note that the abundance of Si, which is the least depleted 
of the three elements discussed here, has the steepest gradi- 
ent with < n(H) >. This reflects the same, unexpected result 
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Figure 6. Relative Ti/O abundance as a function of D/H. Re- 
sults from this work are show by stars; filled circles are data 
points taken from Prochaska et al. (2005). The shaded region 
shows the constraints on the primordial value of D /H determined 
from quasar absorption line systems. For two elements with a sim- 
ilar nucleosynthetic origin, but different refractory properties, we 
expect a correlation with D/H if depletion onto dust is the cause 
of scatter in deuterium abundances. 

presented in Figure |31 i.e. that the depletion fractions seem 
to be inverted relative to our expectations. This may indi- 
cate that the normal circumstances which govern the corre- 
lation between depletion and density in the ISM at large is 
not at work in these local (mostly d < 500 pc) sightlines. 
However, we should also stress that we have many fewer 
sightlines than the Jenkins sample, so that small number 
statistics could contribute to the absence of a tight corre- 
lation. Also noteworthy is the lack of correlation, and very 
large scatter at any given density, of [D/hQ. We have only 
plotted in Figure[5]those stars in our Ti II sample. The same 
plot, but with a much larger sample extending to lower val- 
ues of < n(H) >, has been recently presented by Oliveira 
et al. (2006) who note that the scatter in D/H increases 
dramatically when < 7i(H) > > 0.1, the regime in which 
most of our points lie. If lower D/H values do result from 
differential depletion, it appears from our Figure [5] that this 
depletion does not correlate with < n(H) >. Oliveira et al. 
(2006) have attempted to fit an abundance/depletion model 
to their larger dataset and claim an anti-correlation between 
D/H and < n(H) >. However, the reduced for this fit is 
large (5.8), the scatter in D/H for a given < n(H) > is ~ 
1 dex (similar to what is seen in our data) and the anti- 
correlation driven by just 2 data points. Indeed, some of 
the highest D/H values occur at the highest < n(H) >. We 
therefore conclude that there is no convincing evidence for 
an anti-correlation between < n(H) > and D/H as one might 
expect from dust depletion. 

Although the acceptance of dust depletion as the cause 
of the local variations of deuterium is gaining momentum, 
our results indicate that this interpretation is not clear 
cut. The same conclusion has also been recently reached by 
Steigman et al. (2007) who consider the evidence for dust 

* Here, [D/H] represents the abundance relative to the primordial 
value, which we take to be 26 ppm. 



depletion as 'ambiguous'. However, the above discussion as- 
sumes a very simple approach to the behaviour of dust in the 
local ISM. There are several caveats that should therefore 
be applied, such as the many flavours of dust that exist with 
different chemical compositions and the fact that dust pro- 
duction and destruction is likely to be a dynamic process. 
Moreover, current models of deuterium depletion involve D 
incorporation into deuterated PAHs (Draine 2004, 2006), 
different again from the grains which adsorb most refractory 
elements. Combined with the sightline averaged quantities 
that we are necessarily dealing with, some scatter in abun- 
dance relations is perhaps not surprising and we return to 
this issue in section 14.31 Shocks may also disrupt 'normal' 
depletion patterns. For example, weak shocks may be suf- 
ficient to release D from the dust phase (since deuterated 
PAHs are relatively fragile), whereas much stronger shocks 
may be required to ablate Ti grain material. These issues 
highlight the complexity of the ISM and abundance ratios 
and the caution that should be exercised in the interpreta- 
tion that they are driven by dust. 

4.2 Infall 

An alternative scenario to the depletion of deuterium onto 
grains is the infall of unenriched gas from the halo which 
dilutes the abundance of chemical elements. Knauth, Meyer 
& Lauroesch (2006) have recently suggested local infall to 
explain the difference in N/0 abundances at d > 500 pc, 
compared with closer sightlines. Romano et al. (2006) have 
also combined infall with astration to result in the relatively 
modest destruction (/ = (D/H) prim/ (D/H)lism < 1.8) of 
deuterium from its primordial value. In most simple infall 
scenarios, the infalling gas has a D/H abundance that is 
close to the primordial value (i.e. higher than typical ISM 
values), whilst the abundances of other elements are lower 
than in the ambient disk ISM. This would result in an anti- 
correlation between heavy element abundance and D/H (e.g. 
Steigman et al. 2007), rather than the positive correlations 
depicted in Figure |3] Infall of low metallicity material could 
only explain the positive correlations seen in Figure 2] if the 
infalling low metallicity gas is also poor in deuterium. Al- 
though unlikely, we briefly explore this possibility here. 

One might discriminate between the two scenarios of 
dust depletion versus various infall scenarios by examining 
relative abundances of different heavy elements. If infall of 
unenriched (or low metallicity, but with a solar abundance 
pattern) gas is responsible for the spatial variations in D /H, 
the relative abundances of two elements with similar nu- 
cleosynthetic origins should not correlate with D/H. Con- 
versely, if the same two elements have different refractory 
properties, they should vary with D/H if depletion is the 
cause of variations in the deuterium abundance. In Figure 
[6] we compare the abundances of Ti and O, which are both 
Q elements whose enrichment comes mainly through Type 
II supernovae of massive stars. Applying the above logic, 
one would expect a flat distribution of Ti/O vs. D/H in 
the infall scenario, but a positive correlation with D/H for 
dust depletion. The Spearman rank correlation coefficient is 
—0.1, which indicates that there is no significant correlation 
between Ti/O and D/H, and the negative value actually in- 
dicates the data tend to an anti-correlation. Unfortunately, 
this test is not conclusive, since the very mild dependence 
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Figure 5. Abundances versus total hydrogen volume density, n(H) for a sample of local Galactic sightlines from Linsky et al. (2006) and 
this paper. Whereas Ti, Si and Fe are plotted relative to solar, the deuterium abundances are plotted relative to the primordial value, 
which we take to be 26 ppm. The solid lines show a fit to [X/H] vs. < n(H) > determined from 144 stellar sightlines (Jenkins 2004). The 
dashed lines are the least squares fits to the data presented and have equations: [Ti/H]= — 0.3K «(H) > —1.98, [Si/H]= — 1.28< fi(H) > 
-1.86, [Fo/H]= -0.39< n(H) > -1.75. 



of Ti/H on D/H (Figure |4| spans a factor of three less (0.5 
dex) than the scatter in the Ti/O ratios at a given D/H. 



4.3 Scatter 



We have already argued that a 'standard' infall scenario 
with low metallicity, enhanced D /H gas would cause an anti- 
correlation between [X/H] and D/H. Although D-poor infall 
could qualitatively explain the positive correlations in Fig- 
ure |4l again the gradients provide strong evidence against 
this possibility. If infalling gas has a low metallicity, but a 
solar abundance pattern, then in a plot of [X/H] vs. D/H, 
every element should have the same gradient (ignoring all 
other effects such as depletion, ionization etc.). This is sim- 
ply due to the fact that each element is diluted by the same 
amount and with a logarithmic ordinate, the Ay is a mea- 
sure of the dilution factor. The very different slopes seen 
in Figure |4] therefore argue against a D-poor infall model, 
although we can not rule out more contrived models with 
highly non-solar abundances. However, combined with other 
evidence against the infall model, such as the relative con- 
stancy of certain elements such as O and Kr (Oliveira et al. 
2005; Cartledge et al. 2004) and the sharp transition from 
approximately constant D/H inside the local bubble to a 
large scatter at distances beyond 100 pc, which would re- 
quire a coincidentally very localized region of infall (e.g. de 
Avillez & Mac Low 2002) , we conclude that infall is unlikely 
to be the main cause of local D/H variations. This is not to 
say that the local ISM is a closed box, in fact Romano et al. 
(2006) require some infall to explain the modest astration 
between the primordial D abundance and the local values. 



In closing this section, we note that for a given D/H, there 
is a wide range of observed [Ti/H] abundances, spanning 
up to a factor of three (Figure [J). The scatter means that 
some D/H~10 ppm sightlines have [Ti/H] as high as the 
highest D/H sightlines, and vice versa. Linsky et al. (2006) 
discussed a number of factors that might contribute scat- 
ter in the relationship between element depletion and D/H 
such as line saturation of DI transitions, ionization correc- 
tions and N(HI) determination. The former of these may 
affect our results, although when saturation is obvious, we 
remove the sight line from our analysis (e.g. HD 206773). 
However, we expect the latter two concerns to be alleviated 
by our choice of high N(HI) sight lines. The more obvious 
damping wings should facilitate the determination of N(HI) 
that can be more challenging for log N(HI) < 19 and the gas 
will be approaching complete self-shielding at these high col- 
umn densities. In short, the observed scatter in Ti/H is most 
likely to represent true variations in the gas-phase abun- 
dance, not systematic effects. Such intrinsic scatter may be 
expected due to the dynamic nature of the ISM as dust is 
formed and ablated and the gas phase responds to photoion- 
ization and even incorporation into molecules. Steigman et 
al. (2007) have also argued that some amount of infalling 
gas that is not well mixed can not be ruled out. This will 
also contribute to intrinsic scatter in abundances. 
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5 CONCLUSIONS 

We present 12 new measurements of Ti/H for local Galactic 
sightlines for which D/H measurements exist, almost tripling 
the sample from Prochaska et al. (2005). We investigate el- 
emental abundance trends in the context of depletion of D 
onto grains as the cause for the large variation soon in D/H 
in the local Galactic disk (e.g. Wood ct al. 2004 and refer- 
ences therein). Our data support the tentative correlation 
between Ti/H and D/H found by Prochaska et al. (2005), 
showing a correlation at 97% significance. However, it is not 
clear that this correlation can bo explained by a simple de- 
pletion model. This conclusion is based on 1) the much shal- 
lower dependence of Ti/H on D/H than Si or Fe, despite 
the fact that the former is much more refractory and 2) the 
lack of a steep correlation between [Ti/fi] and < n(H) >. 
If the variation in Ti abundances of our sample sightlines 
were dominated by varying dust depletion, we would expect 
a very steep correlation with < n(H) >, as has been pre- 
viously found (Gondhalekar 1985; Jenkins 2004). Moreover, 
the correlation of Si/H versus D/H is steeper than between 
Fe/H and D/H, contrary to what is expected if these rela- 
tions trace depletion, since Fe is more refractory than Si. 
This may mean that whatever effect is causing the flatten- 
ing of the Ti depletion correlations that arc normally seen in 
the ISM at large, may also affect Fe and Si, although prob- 
ably to a lesser degree. However, the complex and dynamic 
nature of the ISM and dust depletion no doubt complicates 
the above simplistic expectations. Nonetheless, the shallow 
abundance correlations exhibited by titanium are somewhat 
surprising given previous ISM studies, and it is unclear why 
the depletion correlations have broken down in this sample. 
It may be due to superposition of the local bubble compo- 
nents whose physical environment could be quite different 
from the predominantly neutral disk ISM. However, only 
using stars at distances > 400 pc does not change our re- 
sults and many stars in the Jenkins (2004) sample are at 
distances less than 500 pc. Finally, we have also argued that 
the gradients of [Si/H], [Fe/H] and [Ti/H] vs. D/H do not 
support a simple model of infall of metal-poor gas with solar 
abundances as the dominant cause of local D variations. 

In summary, whilst the depletion of D onto grains re- 
mains a plausible explanation for the variation in D/H in 
the local ISM, our results show that this interpretation is 
far from clear cut. 
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